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Abstract 

This  paper  presents  the  implementation  of  a  technique  to  minimize  the  input  ripple  current  in  three-phase  voltage- source  pulse  width 
modulated  (PWM)  inverters  when  supplying  unbalanced/nonlinear  loads.  This  is  achieved  without  increasing  the  KVA  ratings  of  the  dc-link 
filter  components.  The  results  presented  are  particularly  important  for  designers  of  fuel  cell  power  systems.  Inverter  input  ripple  current  has 
been  reported  to  possibly  degrade  fuel  cell  performance  and  reduce  operating  life  if  ripple  currents  are  not  adequately  controlled.  Experimental 
results  from  a  prototype  are  included  to  verify  the  findings  from  the  analysis  and  simulation. 

©  2005  Elsevier  B.Y.  All  rights  reserved. 
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1.  Introduction 

Fuel  cell  systems  are  being  considered  for  utility  power 
applications  or  for  on-site  generation  for  domestic/industrial 
use  [1,2].  One  complication  is  that  fuel  cells  generate  direct 
current  (dc)  while  the  electric  power  system  and  thus  existing 
end  use  equipment  has  been  designed  for  alternating  current 
(ac).  Therefore,  the  fuel  cell  associated  power  electronics 
system  must  include  a  dc  to  ac  converter  (inverter).  Inverter 
input  ripple  current  has  been  reported  to  possibly  degrade 
fuel  cell  performance  and  reduce  operating  life  if  ripple  cur¬ 
rents  are  not  adequately  controlled  [3,4].  Since  the  reactant 
utilization  is  known  to  impact  the  mechanical  nature  of  a 
fuel  cell,  it  is  suggested  in  Ref.  [3]  that  the  varying  reactant 
conditions  surrounding  the  cell  (due  to  ripple  current) 
govern,  at  least  in  part,  the  lifetime  of  the  cells.  In  Ref.  [4], 
it  has  been  shown  experimentally  that  the  ripple  current  can 


*  Corresponding  author. 

E-mail  addresses:  wshireen@uh.edu  (W.  Shireen),  arefeen@ieee.org 
(M.  Arefeen). 

0378-7753/$  -  see  front  matter  ©  2005  Elsevier  B.V.  All  rights  reserved, 
doi:  10. 1016/j. jpowsour. 2005. 06.012 


contribute  to  a  reduction  in  the  fuel  cell  available  output 
power,  cause  internal  losses  and  increase  distortion  of  its 
terminal  voltage.  Passive  filters  (inductor  and/or  capacitor) 
are  typically  specified  to  reduce  ripple  currents  in  the  dc-link 
between  the  fuel  cell  and  the  inverter  to  “perceived”  low 
risk  values.  In  order  to  meet  the  system  requirements,  the 
common  practice  is  to  use  conservative  passive  filters  in  the 
dc-link  that  are  bulky,  expensive  and  inherently  unreliable. 
In  a  recent  paper,  the  authors  have  focused  on  developing  an 
equivalent  circuit  model  of  a  fuel  cell  to  evaluate  the  effects 
of  inverter  input  ripple  current  on  the  performance  of  a  fuel 
stack  [4]. 

Advanced  pulse  width  modulation  (PWM)  techniques  are 
used  to  control  the  semiconductor  switching  devices  in  fuel 
cell  inverters,  in  order  to  meet  the  output  voltage  and  fre¬ 
quency  specifications  [5,6].  Most  advanced  PWM  methods 
used  for  inverter  control  ensure  minimum  harmonics  at  the 
output.  The  inverter  input  current  ripple  and  input  power  fac¬ 
tor  have  also  been  a  point  of  concern  and  various  approaches 
have  been  proposed  so  that  the  total  harmonic  distortion 
(THD)  of  the  inverter  input  current  can  be  kept  within  the 
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desired  limit  in  accordance  to  the  IEEE  519  standards  [7-10] . 
The  approaches  in  Refs.  [10,1 1]  consider  balanced  and  linear 
three-phase  loads.  However,  a  large  majority  of  new  loads  are 
nonlinear  in  nature.  Also,  single-phase  loads  connected  to  a 
three-phase  system  in  most  cases  result  in  unbalanced  load 
conditions. 

In  view  of  this,  this  paper  presents  the  analysis  of  a 
three-phase  PWM  inverter  system  focusing  on  assessment 
of  the  inverter  input  current  with  balanced,  unbalanced  and 
nonlinear  loads.  Analytical  equations  using  the  switching 
function  approach  are  used  to  find  the  proper  state  equations 
to  describe  the  power  conversion  circuit  in  MATLAB  [12]. 
Recent  research  has  shown  that  the  switching  function 
concept  is  a  powerful  tool  in  the  analysis  of  static  power  con¬ 
verters  [13].  In  this  paper,  the  switching  function  approach 
used  to  derive  the  state  equations  is  described  and  simulation 
results  from  the  MATLAB  program  are  presented.  Effects 
of  balanced,  unbalanced  and  nonlinear  loads  on  the  inverter 
input  current  are  presented.  It  is  found  that  unbalanced  linear 
or  nonlinear  loads  at  the  output  contribute  to  the  appearance 
low  frequency  ripple  component  in  the  inverter  input  current. 
These  ripple  currents  must  be  supplied  from  the  dc-link, 
which  necessitates  an  increase  in  the  KVA  ratings  of  the 
dc-link  filter  components.  Based  on  the  analytical  findings, 
a  modified  PWM  method  is  proposed  to  minimize  the  lower 
order  ripple  currents  at  the  inverter  input  due  to  unbalanced 
load  currents  without  increasing  the  KVA  rating  of  the  dc- 
link  filter  capacitor.  The  proposed  method  requires  sensing 
of  the  load  currents  and  balancing  out  the  power  in  the  three 
phases  even  when  the  loads  are  unbalanced  by  modifying 
the  inverter  modulation  in  small  steps.  Experimental  results 
are  included  to  verify  the  findings  from  the  analysis  and 
simulation. 

Reliability  and  cost  of  the  power  electronics  system  are 
important  considerations  for  the  commercialization  of  fuel 
cell  power  sources.  Inverter  input  ripple  current  has  been 


reported  to  possibly  degrade  fuel  cell  performance  and  reduce 
operating  life.  The  analytical  method  and  proposed  ripple 
minimization  technique  will  enable  fuel  cell  system  designers 
to  incorporate  more  reliable  power  electronics  in  their  design 
by  taking  into  account  inverter  ripple  current  due  to  modern 
nonlinear/unbalanced  loads. 


2.  Analysis  of  three-phase  PWM  inverters — the 
switching  function  approach 

Fig.  1  shows  a  typical  three-phase  PWM  inverter  sys¬ 
tem.  The  dc-input  for  the  three-phase  inverter  is  obtained 
by  rectifying  a  60  Hz  ac  source  by  means  of  a  diode  rectifier. 
The  dc-link  filter  capacitor  are  so  selected  that  the  voltage 
Vdc  at  the  inverter  input  is  constant  and  ripple-free.  The 
semiconductor  switching  devices  (S  i— S6)  of  the  inverter  are 
controlled  by  PWM  signals  to  obtain  three-phase  near  sinu¬ 
soidal  ac  voltages  of  the  desired  magnitude  and  frequency  at 
the  inverter  output. 

In  order  to  analyze  the  PWM  inverter  system,  it  is  impor¬ 
tant  to  compute  a  dependent  variable  (input  current  and  output 
voltage)  in  terms  of  an  independent  variable  (input  voltage 
and  output  current).  The  switching  function  approach  is  used 
to  develop  a  functional  model  of  the  inverter  that  can  be  easily 
simulated  in  MATLAB. 

For  a  three-phase  inverter,  the  transfer  function  (7)  is  com¬ 
posed  of  three  independent  switching  functions  and  is  given 

by, 

r  =  [sfl_isfo_isc_i]  (i) 

The  switching  functions  are  Fourier  series  representation 
of  the  switching  sequence  used  for  PWM  control  of  the 
inverter  switching  devices.  Mathematical  representation  of 


Fig.  1.  Three-phase  PWM  inverter  system. 
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Fig.  2.  SPWM  generation. 


Fig.  3.  Switching  function  Sa-\. 


the  switching  functions  are  given  by, 


oo 

Sa-i(a)t )  =  22An  Sb-i(a>t) 

n— 1 


Sc-i(cot) 


The  above  switching  functions  can  be  generated  by  using 
the  sinusoidal  pulse  width  modulation  (SPWM)  technique 
as  shown  in  Fig.  2.  The  ‘saw  tooth’  function  in  MATLAB 
was  used  to  generate  the  carrier  signal  of  the  desired  fre¬ 
quency.  The  carrier  signal  frequency  used  for  the  simulation 
was  1980  Hz.  The  switching  function  Sa- 1  for  phase  V  of 
the  inverter  was  generated  by  comparing  the  carrier  signal 
with  modulating  signal  Mod_a.  Similarly,  switching  func¬ 
tions  Sb- 1  and  Sc-\  can  be  obtained  by  comparing  carrier 
signal  with  modulating  signals  Mod  J?  and  Mod_c  for  phases 
‘b’  and  ‘c’,  respectively.  Fig.  3  shows  the  switching  function 
Sa- 1 .  The  frequency  spectrum  of  any  periodic  signal  can  be 
obtained  by  using  the  fast  Fourier  transform  (FFT)  function 
in  MATLAB.  It  can  be  seen  from  Fig.  4  that  the  switching 
function  consists  of  the  fundamental  60  Hz  component  and 
higher  order  switching  frequency  harmonics.  Lower  order 
harmonics  up  to  the  33rd  (1980  Hz)  has  been  eliminated  by 
the  SPWM  strategy. 

The  switching  functions  Sa- 1,  Sb- 1  and  Sc- 1  are  used  to 
calculate  the  inverter  output  voltages.  The  voltages  vao ,  Vbo 
and  vco  can  be  obtained  as, 


[Vao  Vbo  Fco] 


where,  Vdc  is  the  inverter  input  voltage.  Also, 

(' vao  +  vbo  +  Vco) 

Vno  —  ~ 


Fig.  4.  Frequency  spectrum  of  Sa-\. 


Hence,  assuming  star  connected  loads  at  the  output  of  the 
inverter,  the  inverter  line  to  neutral  output  voltages  are  given 

by, 


Van  —  Va0  Vn0,  Vbn  —  Vbo  Vn0,  Vcn  —  Vco  Vbo- 


3.  Analysis  of  input  current  ripple  using  MATLAB 

For  a  given  load  condition,  the  inverter  input  current  can 
be  calculated  by  using  another  set  of  three-phase  switching 
functions  Sa- 2,  Sb- 2  and  Sc~ 2-  The  switching  function  Sa- 2  is 
similar  to  Sa- 1 ,  the  only  difference  is  that  Sa- 1  varies  from  —  1 
to  +1,  whereas  Sa- 2  varies  between  0  and  +1.  The  frequency 
spectrum  for  Sa- 2  will  be  similar  to  5^-1  (in  Eq.  (2))  with  an 
additional  dc  component.  The  inverter  input  current  U  can  be 
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calculated  as, 

Ii(cot)  =  Ia(cot)Sa- 2  +  Ib(cot)Sb- 2  +  Ic(cot)Sc- 2  (6) 

where,  Ia(tot),  h(tot)  and  Ic(tot)  are  the  three-phase  load  cur¬ 
rents. 


3.1.  Inverter  input  current  with  balanced  linear  load 


If  balanced  three-phase  linear  loads  of  impedance  Z(co) 
are  connected  at  the  inverter  output,  then  the  inverter  output 
line  currents  can  be  computed  as, 


Ia(cot)  = 


Van  (.tot ) 

Z((o)  ’ 


lb(cot)  —  I  a 


Ic((Ot)  —  I a 


The  mathematical  operation  in  Eqs.  (6)  and  (7)  were  per¬ 
formed  in  MATLAB  to  obtain  the  inverter  input  current  for 
balanced  loads  and  its  corresponding  frequency  spectrum 
is  shown  in  Fig.  5.  From  Fig.  5,  it  can  be  seen  that  the 
inverter  input  current  (//)  with  balanced  loads  consists  of  a  dc 
component  and  higher  order  switching  frequency  harmonics 
(1980  Hz  and  above). 
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Fig.  6.  FFT  of  7/  with  unbalanced  load. 


The  same  approach  used  for  balanced  loads  can  be  used  to 
plot  the  inverter  input  current.  Fig.  6  illustrates  the  FFT  of  the 
inverter  input  current  with  unbalanced  loads.  It  can  be  seen 
that  in  addition  to  the  dc  component,  the  inverter  input  cur¬ 
rent  consists  of  120  Hz  component  of  significant  magnitude. 
Hence,  unbalanced  linear  loads  at  the  inverter  output  generate 
abnormal  lower  order  harmonics  at  the  inverter  input  which 
are  not  present  in  the  case  of  balanced  loads. 

3.3.  Inverter  input  current  with  balanced  nonlinear  load 


3.2.  Inverter  input  current  with  unbalanced  linear  load 


If  three-phase  unbalanced  loads  Za(co),  Zb(co)  and  Zc(co) 
are  connected  at  the  inverter  output,  using  the  voltages  from 
Eq.  (5),  the  three-phase  unbalanced  currents  can  be  calculated 
as, 


Ia(cot)  = 


h(tot)  = 


Ic(cot)  = 


O an  (tot) 

Za 

Vbn(tOt) 


Zb 

Vcn  (tOt) 


Frequency,  in  Hz 


Fig.  5.  FFT  of  f  with  balanced  load. 


When  loads  connected  at  the  output  of  the  inverter  are 
nonlinear  in  nature,  the  load  currents  consists  of  harmonics  in 
addition  to  the  fundamental  frequency  component.  Assum¬ 
ing  the  following  balanced  but  nonlinear  load  currents  are 
assumed  for  inverter  input  current  calculation. 


Ia(cot)  =  In  sin (ncot 

n— 1,3, 5, 7, 9 


Ib(tOt)  —  I a 
3c  (tot)  =  la 


Using  Eq.  (6),  the  inverter  input  current  was  calculated.  The 
FFT  of  the  input  current  showed  the  appearance  of  lower 
harmonics  (6th,  12th,  etc.)  in  addition  to  the  dc  component  in 
the  inverter  input  current  spectrum.  But,  these  harmonics  are 
very  small  in  magnitude.  Hence,  balanced  nonlinear  loads  at 
the  inverter  output  do  not  produce  any  significant  lower  order 
harmonics  at  the  inverter  input. 


3.4.  Inverter  input  current  with  unbalanced  nonlinear 
load 


Fig.  7  shows  the  FFT  of  7/  with  unbalanced  nonlinear 
loads  at  the  inverter  output.  Note  the  appearance  of  all  the 
abnormal  lower  order  harmonics  (2nd,  4th,  5th,  etc.).  Hence, 
unbalanced  nonlinear  loads  at  the  inverter  output  generate 
abnormal  lower  order  harmonics  at  the  inverter  input. 
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Fig.  7.  FFT  of  7/  with  unbalanced  nonlinear  load. 

4.  Minimization  of  input  ripple  current 


when  n=  1,  but  we  also  get  a  second  harmonic  component 
of  magnitude  IaA\ .  Similar  components  are  obtained  from 
the  products  Ib(cot)Sb- 2  and  Ic(cot)Sc- 2.  When  the  loads  are 
balanced,  Ia=Ib=Ic  and  0ioad-a  =  0ioad-b  =  0ioad-c.  Hence, 
the  second  harmonic  components  are  also  balanced  and  the 
sum  of  the  three  components  adds  up  to  zero.  When  the  load 
currents  are  unbalanced,  the  second  harmonic  current  com¬ 
ponent  obtained  from  the  product  terms  in  Eq.  (6)  do  not  add 
up  to  zero  and  that  result  in  the  appearance  of  the  second  har¬ 
monic  in  the  inverter  input  current.  Based  on  this  analysis, 
the  proposed  technique  involves  phase  correction  and  mag¬ 
nitude  correction  of  the  modulating  signals  used  for  inverter 
control.  This  would  ensure  that  sum  of  the  second  harmonic 
components  resulting  from  the  three  product  terms  in  Eq.  (6) 
add  up  to  zero  regardless  of  the  unbalance  in  the  load. 

4.1.  Phase  angle  correction 


It  is  evident  from  the  analysis  presented  in  the  earlier  sec¬ 
tions  that  the  inverter  input  current  contains  low  frequency 
ripple  components  if  unbalanced/nonlinear  loads  are  con¬ 
nected  at  the  inverter  output.  As  discussed  in  Section  1, 
these  abnormal  harmonics  must  be  supplied  from  the  de¬ 
link,  which  necessitates  an  increase  in  the  KVA  ratings  of 
the  dc-link  filter  components.  In  this  section,  a  ripple  mini¬ 
mization  technique  is  proposed  which  does  not  necessitate  the 
increase  in  dc-link  filter  values  even  with  unbalanced  loads 
at  the  inverter  output. 

The  three-phase  switching  functions  Sa- 2,  Sb- 2  and  Sc- 2 
mentioned  in  Section  3  can  be  mathematically  represented 
as, 


00 

Sa-i(0)t)  =  A0  +  EA  sin (n(cot  +  9a)), 

n  =  1 


00 

Sb-2(o)t)  =  Ao  +  'y  ] An  sin 

n  =  1 
00 

Sc-2((Dt )  =  A0  +  EA  sin 

n= 1 


(10) 


where,  0a ,  0b  and  0C  are  phase  shifts  in  the  three-phase 
modulating  signals.  Under  normal  balanced  conditions, 

oa=ob=oc= 0. 

Let  the  three-phase  output  currents  of  the  inverter  be, 


—  Ia  sin(<Z>7  $load— a)-> 

f  271 

Ib((i)t )  =  Ib  sin  (  COt  $load— b 


(ID 


2tx 

Ic((Pt)  —  Ic  sin  I  cot  A  ~z~  $load- 


where,  0\oad-a ,  0\oad-b  and  0\oad-c  represent  the  impedance 
angles  associated  with  the  three-phase  loads. 

While  calculating  the  input  current  using  Eq.  (6),  the  prod¬ 
uct  of  Ia(ojt)  and  Sa- 2  results  in  harmonic  components  at 
frequencies  (n  —  1  )co  and  (n  +  l)co.  We  obtain  a  dc  component 


Substituting  Eqs.  (10)  and  (11)  in  Eq.  (6)  gives  the  input 
current  expression.  Using  the  input  current  expression,  the 
phase  angle  correction  that  needs  to  be  introduced  in  the 
modulating  signals  of  each  phase  in  order  to  make  the  sec¬ 
ond  harmonic  components  in  the  input  current  expression 
balanced  is  given  by, 


$load —b 


(12) 


4.2.  Magnitude  correction 

The  inverter  modulation  needs  to  be  modified  to  maintain 
the  equality, 


MaIa  =  MbIb  =  MCIC  (13) 

where,  Ma,  Mb  and  Mc  are  the  magnitudes  of  the  modulat¬ 
ing  signals  Mod_a,  ModJ?  and  Mod_c  shown  in  Fig.  2.  The 
following  is  the  step-by-step  process  to  implement  the  phase 
and  magnitude  correction  technique: 

1.  Initially,  modulating  signals  are  equal  in  magnitude  and 
are  displaced  120°  from  each  other. 

2.  Load  currents  are  measured,  if  they  are  unbalanced,  phase 
shifts  equal  to  half  of  the  impedance  angles  are  introduced 
to  the  modulating  signals  of  the  corresponding  phase. 

3.  Load  currents  are  measured  and  checked  for  the  equality 
in  Eq.  (13). 

4.  Ma,  Mb  and  Mc  are  changed  in  small  steps,  keeping  their 
average  value  same  as  that  in  step  [1]. 

5.  Steps  [3]  and  [4]  are  repeated  until  the  desired  equality  in 
step  [3]  is  reached. 

Fig.  8  shows  the  simulated  inverter  input  current  with  the 
above  technique  in  effect.  Comparing  with  Fig.  6,  it  can  be 
seen  that  the  2nd  harmonic  component  has  been  reduced  to 
nearly  zero. 
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Fig.  8.  FFT  of  f  with  unbalanced  load  (with  correction). 


5.  Experimental  results 


Tek  Run:  2.50kS/s  Sample  IliPM 


A:1. 075kHz 
@: 103  Hz 

Cl  Mean 
1.816  V 


Chi  2.00  VQ 

400mV 
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80  mV  9  Apr  2004 
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A  laboratory  prototype  of  a  three-phase  inverter  as  shown 
in  Fig.  1  was  implemented  to  verify  the  findings  from  the  anal¬ 
ysis  and  simulation.  The  dc-input  voltage  was  set  at  50  V.  A 
DSP  controller  from  Texas  Instruments  C2000  family  was 
used  to  control  the  system.  The  DSP  controller  provided  on- 
chip  PWM  modules  with  multiple  PWM  outputs.  The  on-chip 
10-bit  analog-to-digital  converter  was  utilized  for  data  acqui¬ 
sition  to  obtain  load  current  information.  The  16-bit  fixed 
point  DSP  core  was  the  computational  engine  for  the  imple¬ 
mentation  of  the  proposed  technique.  The  PWM  switching 
frequency  was  set  at  20  kHz  with  a  fundamental  frequency  of 
60  Hz.  Fig.  9  shows  the  experimental  inverter  input  current 
with  balanced  RL  loads.  Lamp  loads  were  used  as  the  resis¬ 
tive  load.  The  FFT  of  the  current  in  Fig.  9  shows  the  presence 
of  the  dc  component  and  higher  order  normal  switching  fre¬ 
quency  harmonics,  which  agrees  with  the  simulation  results. 
Then,  the  inverter  prototype  was  tested  with  unbalanced  RL 
loads.  This  was  achieved  by  connecting  different  number  of 
lamp  loads  in  each  phase  while  keeping  the  inductor  in  each 


Fig.  10.  Experimental  /;  and  FFT  with  unbalanced  loads  (without  correc¬ 
tion). 


Mathl 
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Fig.  11.  Experimental  /;  and  FFT  with  unbalanced  loads  (with  correction). 
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phase  constant.  Fig.  10  shows  the  experimental  inverter  input 
current  with  unbalanced  loads.  It  can  be  seen  that  the  input 
current  now  consists  of  a  2nd  harmonic  (120  Hz)  component 
of  significant  magnitude  (about  30%  of  the  dc  component). 
Fig.  1 1  shows  the  inverter  input  current  with  the  same  unbal¬ 
anced  loads  (as  in  Fig.  10)  but  with  the  proposed  correction 
technique  in  effect.  It  can  be  seen  that  the  lower  order  2nd 
harmonic  (120  Hz)  component  has  been  significantly  reduced 
without  affecting  the  dc  component  of  the  current. 


6.  Conclusions 

A  ripple  minimization  technique  to  minimize  the  input 
ripple  current  in  three-phase  voltage- source  pulse  width 
modulated  inverters  when  supplying  unbalanced/nonlinear 
loads  is  presented.  The  proposed  technique  can  contribute 
in  improving  the  reliability  and  reducing  the  cost  of  the 


Fig.  9.  FFT  of  f  with  balanced  load. 
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power  electronics  system  in  the  commercialization  of  fuel 
cell  power  sources.  The  proposed  technique  is  verified  in  an 
actual  setup  and  experimental  results  are  included  to  verify 
the  findings  from  analysis  and  simulation. 
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